A Hydrodynamic Description of Heavy Ion Collisions at the SPS and RHIC 
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A hydrodynamic + cascade model of relativistic heavy ion colhsions is presented and com- 
pared to available hadronic data from the SPS to RHIC. The model consistently reproduces the 
radial and elliptic flow data for different particles, collision energies, and impact parameters. 
Three ingredients are essential to the success: (a) a reasonable EOS exhibiting the hard and soft 
features of the QCD phase transition, (b) thermal hadronization at the phase boundary, and 
(c) subsequent hadronic rescattering. Some features of the RHIC data are readily explained: (i) 
the observed elliptic flow and its dependence on pt and mass, (ii) the anomalous p/n~ ratio for 
Pt ~ 2.0 GeV, (iii) the difference in the slope parameters measured by the STAR and PHENIX 
collaborations, and (iv) the respectively strong and weak impact parameter dependence of the 
p and (j) slope parameters. For an EOS without the hard and soft features of the QCD phase 
transition, the broad consistency with the data is lost. 



I. INTRODUCTION 

A. Reaching the Macroscopic Limit 

Excited nuclear matter has been created by colliding 
Pb and Au ions at the SPS {^/s = 17 GeV per nucleon) 
and RHIC {y^ = 130 GeV per nucleon) accelerators 
0, H, 1^ . For infinitely large nuclei, the excited matter can 
be characterized by macroscopic quantities - pressure, 
temperature, viscosity, etc. Lattice QCD simulations in- 
dicate that for temperatures larger than Tc « 160 MeV, 
confined nuclear matter morphs into a phase of decon- 
fined quarks and gluons - the Quark Gluon Plasma 
(QGP) I, |. The possibihty of observing the QGP in 
real nuclei has motivated the heavy ion experimental pro- 
gram. 

If the system is macroscopic then thermodynamics de- 
scribes the static properties of the matter and hydrody- 
namics describes the dynamic properties of the matter. 
In fact, the observed particle ratios are remarkably close 
to the particle ratios in an ideal gas of hadrons at a tem- 
perature, T « 165 MeV §, 0, |, |. This suggests that 
the system evolved from a state close to thermal equi- 
librium at the phase transition boundary. However, the 
same thermal description reproduces the hadron ratios in 
proton-proton and e+e^ collisions, where the system size 
is small and equilibration seems impossible. The success 
of thermodynamics seems to reflect phase space rather 
than the equilibration of macroscopic system. 

Therefore, a thermodynamic (static) description, di- 
vorced from a hydrodynamic (dynamic) description, can 
not unambiguously signal a macroscopic state. It is im- 
portant that elementary proton-proton and e'^e^ colli- 
sions do not exhibit hydrodynamic behavior. An analy- 
sis of hadronic spectra shows little sign of the trans- 
verse expansion predicted by hydrodynamics. Thus, the 
excited systems produced in these elementary collisions 



are not macroscopic. 

In contrast, experiments with heavy ions do show evi- 
dence for a hydrodynamic expansion. Momentum corre- 
lations, colloquially known as flow, are observed at the 
SPS and RHIC. In PbPb collisions, the particles emerge 
from the collision with a collective transverse velocity 
of approximately 1/2 c |jll[. This radial flow is firmly 
established from a combined analysis of particle spec- 
tra, HBT correlations, and deuteron coalescence |jl^. In 
non-central collisions, the particles emerge with an ellip- 
tic flow (see for example[|[ ^). Elliptic flow is quantified 
by V2, the asymmetry of the angular distribution 



V2 = (cos(2(/))) 



(1) 



where (f) is measured around the beam axis with respect 
to the impact parameter. Radial and elliptic flow data 
are measured as a function of transverse momenta, par- 
ticle type, impact parameter, and collision energy. This 
wealth of momentum correlations severely constrains vi- 
able models of the heavy ion collision. 

Several microscopic models have been used to explain 
the available heavy-ion data. The first is a dilute par- 
ton model, which is quantified with the HIJING event 
generator p3[ |. Dilute parton models are based upon the 
extrapolation of perturbative QCD from high px down to 
a scale of ~ 1 GeV. For central AuAu collisions at RHIC, 
HIJING predicted a mini-jet multiphcity of 200, 
which is insufficient to generate the strong hydrodynamic 
response observed at RHIC The second is a string 
model, which is quantified with the UrQMD event gener- 
ator. In string-models non-interacting strings decay into 
hadrons which subsequently interact. Due to the small 
transverse pressure at early times, UrQMD predicted a 
decrease in elliptic flow from the SPS to RHIC A 
« 50% increase was observed. In contrast to these mi- 
croscopic models, hydrodynamic calculations at the SPS 
and RHIC, give a good description of the observed radial 
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and elliptic flows |T|, |T|, |20[ but offer no 
insight into the microscopic mechanism of equilibration. 

Accepting the macroscopic approach and its limita- 
tions, the phase transition to the QGP influences both 
the radial and elliptic flows. Lattice simulations in- 
dicate 1^ that over a wide range of energy densities 
e = 0.5 — 1.3GeV/fm , the temperature and pressure 
are nearly constant and the speed of sound is approx- 



imately zero, 



de 



0. Because the speed of sound is 



small in this range, the pressure can not effectively ac- 
celerate the matter However, when the initial 
energy density is well above the transition region, the 
matter enters the hard QGP phase. The speed of sound 
approaches -^/l/S and the pressure drives collective mo- 
tion. At a time of ~ 1 fm/c, the energy density at the 
SPS and RHIC are very approximately 4 and 8 GeV/fm^ 
p^ , Based on these experimental estimates, the 
hard QGP phase is expected to live significantly longer 
at RHIC than at the SPS. The final radial and elliptic 
flows of the produced particles should reflect this differ- 
ence 

Certainly, hydrodynamics is not applicable when the 
particles decouple from the collision and this "freeze- 
out" must be modeled in order to compare the observed 
hadron spectra to hydrodynamic calculations. Usually, a 
naive freezeout prescription is taken: A "freezeout Tem- 
perature" T/, is specified; thermal and chemical equilib- 
rium are assumed for T < Tf, the spectrum of parti- 
cles passing through the Tj isotherm is calculated; Tf 
is finally adjusted to match the single particle spectrum 
of pions and nucleons. Of course, this prescription is 
unrealistic and takes away from the predictive power of 
hydrodynamics. Nevertheless, the approach successfully 
describes many radial and elliptic observables from the 
SPS to RHIC 0, |l^, 1^. 

However, the naive freezeout prescription fails in a 
number of respects. First, on the time scale of the 
collision ^ 10 fm/c, hadronic reactions do not alter 
the hadron composition and chemical equilibrium is not 
maintained (see for example |3^, |3^, Q). Therefore in 
the late hadronic stages, chemical freezeout must be mod- 
eled in order to describe the resonance contribution and 
the particle ratios. Second, different particle types freeze- 
out at different times and with different transverse veloc- 
ities. With a universal freezeout temperature, the trans- 
verse flow of the strange particles A, S, 17, is never re- 
produced |34|. Third, at the SPS and RHIC, integrated 
elliptic flow is a strong function of the freezeout tem- 
perature. When the universal freezeout temperature is 
adjusted to match the nucleon spectrum, the integrated 
pion elliptic flow is too large |^5|. In reality, pions and 
nucleons freezeout at different times and temperatures. 

To model freezeout, Bass and Dumitru |Q replaced 
the hadronic phase of the hydrodynamics with a hadronic 
transport model, UrQMD. In this approach, the switch 
from hydro to cascade is made at a switching temper- 



the attendant theoretical problems are ignored. The ap- 
proach worked. Chemical freezeout was incorporated 
into a comprehensive dynamical picture. The ffow of 
the multi-strange baryons was reproduced. When a simi- 
lar hydro-f-cascade model |2l[ was applied to non-central 
collisions, elliptic flow was also reproduced at the 20% 
level. Furthermore, with the Tf indeterminacy removed, 
these "simple" boost invariant hydro-f cascade models 
were rather predictive - only ^ and the p/p ratio have 



to be specified. 



B. Brief Summary 



In this work, we compare the hydro-|-cascade model 
of to the fiow systematics at the SPS and RHIC. 
The model uses hydrodynamics to model the initial stage 
of the collision, and the hadronic cascade, Relativistic 
Quantum Molecular Dynamics (RQMD v2.4), to model 
the final stages of the collision . 

The Equation of State (EOS) is varied systematically 
and results are compared to the whole body of flow data 
from the SPS to RHIC. A family of EOS, labeled by 
the value of the Latent Heat (LH) is constructed; LH4, 
LH8, LH16... denote increasingly soft EOS with latent 
heats 0.4, 0.8, 1.6... GeV/fm'^ respectively. As a limiting 
case, the latent heat is made very large forming LHoo. A 
Resonance Gas (RG) EOS is also studied. For an EOS 
with both the hard and soft features of the QCD phase 
transition, the model is broadly consistent with the body 
of available data. The best overall consistency with the 
data is found with LH8. For an EOS with only hard 
(e.g. RG) or only soft features (e.g. LHoo), the broad 
consistency with the data is lost. 

The model consists of three distinct components. The 
flrst component solves the equations of relativistic hy- 
drodynamics in the transverse plane, assuming Bjorken 
scaling |3^. The switching surface, or the isotherm where 
Tswitch = 160 MeV, is found. The sensitivity of the 
model results to TgwUch will be discussed in a separate 
publication where chemical freezeout is also addressed 
||38| . The second component converts the macroscopic 
hydrodynamic variables on the switching surface into 
hadrons according to the Cooper-Frye prescription aug- 
mented with a theta function rejecting backward mov- 
ing particles Finally, the third component, the 
hadronic cascade RQMD sequentially rescatters the gen- 
erated hadrons and models the hadronic freezeout stage 
of the collision. Throughout the analysis the role of 
hadronic rescattering is assessed. In all figures, the "Hy- 
dro-fRQMD" curves incorporate hadronic rescattering 
and resonance decays while the "Hydro Only" curves only 
incorporate resonance decays. 

As outlined in the abstract, several features of the first 
RHIC data are readily explained in the course of this 
analysis. 



ature, Tc, 



litch 



Tc- The spectrum of particles leav- 



ing the surface is taken as the input to the cascade and 
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II. MODEL DESCRIPTION AND THE EOS 



Integrating over the transverse plane, we find that 



A. Hydrodynamics 

Relativistic hydrodynamics is a set of conservation 
laws for the stress tensor (T'"') and for the conserved 
currents (Jf), S^T'"' = and a^Jf = 0. In cquihb- 
rium, T'^'^ and are related to the bulk properties of 
the fluid by the relations, T^"" = {e+p)U''U'' -pg''" and 
Jf = riiU^ Here e is the energy density, p is the 

pressure, tii is the number density of the corresponding 
current, and C/^ — 7(1, v) is the proper velocity of the 
fluid. In strong interactions, the conserved currents are 
isospin (Jj ), strangeness (Jg), and baryon number (Jg)- 
For the hydrodynamic evolution, isospin symmetry is as- 
sumed and the net strangeness is set to zero; therefore 
only the baryon current Jg is considered below. 

The equations of motion may be expressed in terms of 

and r] = jlogd^), which 



the variables r = y/t^ 



are respectively referred to as the Bjorken proper time 
and the spatial rapidity. Boost invariance assumes that 
the solution for any value of 77 may be found by boosting 
the solution at 77 = to a frame moving with veloc- 
ity V = tanh(77) in the negative z-direction. With this 
assumption, the equations of motion become two dimen- 
sional |29, M and are given at 77 = by 



driTT°y) + d^irT^y) + dyirTyy) 



dr{Tj%)+d^{Tjl 



dyirjy^) 



= -p 

= 
= 
= 0. 



(2) 



Integrating over the transverse plane, one finds that net 
baryon number per unit spatial rapidity, / dxdy (tJ^), 
and the transverse momentum per unit rapidity, 
/ dxdy [tT^^) as well as / dxdy (rT'^y), are conserved. 
The energy per unit rapidity, / dxdy (rT^" ) , decreases 
due to the work done per unit time Q by the pressure 
in the longitudinal direction, J dxdy p. 

For an ideal fluid, entropy conservation can be derived 
0, ^f,{S^') = 0. The entropy current is defined as S^' = 
s W^, where s is the entropy density and is the fluid 
4-velocity. For a Bjorken expansion entropy conservation 
becomes 



dx dyrsj = 



dSf, 



drj 



(4) 



is a constant of the motion. This relation is monitored 
to test the accuracy of the solution. 

These equations are solved numerically with a 
Gudunov method 1 42 . Using second order operator split- 
ting Q, a single time step separately updates the x- 
direction, the y-direction, and the loss terms on the r.h.s. 
of Eq. |2[ Different splittings gave only negligibly differ- 
ent results. The simple RHLLE Riemann solver was used 
for the updates in the x and y directions |43, 44]. A sec- 



ond order (in r) Runge-Kutta stepper was used for the 
r.h.s. update. 



B. Initial Conditions 

To model the initial conditions, the entropy and baryon 
distributions at a Bjorken time of tq = 1 fm/c, are made 
proportional to the distribution of participating nucleons 
in the transverse plane. Since entropy and baryon num- 
ber are conserved per unit rapidity, the final yields of 
pions and nucleons are then proportional to the number 
of participants. The initial conditions are similar to sWN 
(entropy per Wounded Nucleon) initial conditions in p^ ] 
and to the initial conditions of . 

For all subsequent discussions, we consider two iden- 
tical (for simplicity) nuclei with atomic number A and 
B, and nucleon distributions pa(^) and psir), collide 
along the z-axis with impact parameter 6, pointing in 
the x-direction from the center of nucleus A, {xa, ys) = 
(—5/2,0), to the center of nucleus B, {xB,yB) = 
(+6/2,0). The nucleon distribution pA is parameterized 



as a Woods-Saxon distribution, p(r) oc 



and 



is normalized to the atomic number A. The parameters 
are 6 = 0.55 fm, Ra = 1.08 A^/^. The Ra used is 4% be- 
low the value used by the STAR collaboration p^^ . The 
number of participating nucleons per unit area, -jiTdy ' 
a position xt — (x, y) in the transverse plane is given by 



dr{TS°) + d^{TS'=) + dy{TSy) = 0. 



(3) 



dNp 
dx dy 



= TA{b/2 + xt) { 1 - 



1 - 



(TNNTB{-b/2 + Xt) 

B 



+ TB{-b/2 + XT) 



1 - 



aNNTAib/2 + Xt) 
A 



(5) 



Here, Ta{xt) = J PA{x^y, z)dz is the thickness of a nu- 
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TABLE I: A summary of the input parameters to the model. 
Cs and C„g are respectively the entropy and baryon number 
per participant per unit rapidity. The values above the double 
line are the input parameters. The values below the double 
hne are derived from the input parameters. The initial energy 
density depends on the EOS and impact parameter. For cen- 
tral collisions and for two EOS spanning the gamut, we quote 
the initial energy density in the center of the collision (eo) and 
the initial energy density averaged over the transverse plane 
with respect to the number of participants ((e)). 



cleus at position (x,y) and a^^N is the inelastic nucleon- 
nucleon cross section. For the sake of comparison, n is 
taken as 33 mb both at the SPS and RHIC. For large A, 
[1 ~ ^ eM-^NNTA{xT)), and often Eq. | 

is re-written in terms of exponents. 

With the number of participants specified, the initial 
entropy and (net) baryon densities at time tq = 1 fm/c, 
are then fixed with two constants Cs and Cng with 



, , Cs dNp 



nB{x,y,To) 



To dx dy 
Cnu dNp 
To dx dy 



(6) 
(7) 



The two dimensionless constants Cg and C„^ are the 
entropy and net baryon number produced per unit spa- 
tial rapidity per participant. At the SPS (see Sect. 
IV A), Cs and C^^ were adjusted to fit the total yield 
of charged particles and the net yield of protons, re- 
spectively. At RHIC, Cs was adjusted to match the 
PHOBOS multiphcity ^ = 555 ± 12(stat) ± 35(syst) 
p5[ . At the, time the p/p ratio was not known and 
sJriB = Cg/CnB 'was estimated from UrQMD simula- 
tions to be w 150. This gives the ratio p/p = 0.45. 
Later, the STAR and PHOBOS collaborations measured 
the ratios, p/p = 0.65 ± .Ol(stat) ± .07(syst) and p/p — 
0.60 ± .04(stat) ± .06(syst) respectively ^ ||. Since 
the measured ratio is close to the ratio initially used, 
and since a full simulation takes several CPU days, the 
UrQMD-based estimate s/ub ~ 150 was used through- 
out this work. This makes the model p yield approxi- 
mately 15% too low and the model proton yield approx- 
imately 15% too high. This correction will be accounted 
for in future works. A summary of the parameters is 
given in Table ^ 

Two quantities, which will be used extensively in the 
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FIG. 1: The derived quantities e and Rrms/RA, defined 
by Eq. ^ and ^ as a function of the number of participants 
relative to the maximum number. The axis on top of the 
graph shows the impact parameter b relative to 2Ra- The 
curves are drawn for PbPb collisions at the SPS, but depend 
only slightly on the colliding system and energy. 



analysis in Sect. [H and Sect. are defined as 

_ (y^-x^) 



{y^ 



(8) 
(9) 



where the average is taken over the initial entropy distri- 
bution of Eq.^. These quantities are plotted as a function 
of the number of participants relative to central collisions 
in Fig. ^ e measures the initial elliptic deformation of 
the overlap region and grows approximately linearly with 
Np. 

For the calculation presented, the entropy and there- 
fore the number of charged particles scales as the number 
of participants. Recently, the experiments have reported 
that the charged particle multiplicit y g rows slightly faster 
than the number of participants |5^. This slight 
growth can be incorporated into hydrodynamics p5[ , but 
instead the experimental dNch/dy is compared directly 
to the model dNch/dy. This makes the model impact 
parameter slightly larger than the impact parameter de- 
termined by the experimental collaborations. 



C. Equation of State 

To solve the equations of motion, we need an Equa- 
tion of State (EOS), or a relation between the pressure 
(p) and the energy and baryon densities (e and n b respec- 
tively). In many previous hydrodynamic calculations, a 
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bag model EOS is used ||l6|, ^ ^ . This has some ad- 
vantages, since the degrees of freedom are explicit in both 
phases. However a typical bag model results in an EOS 
with a large latent heat, LH=1 — 1.5 GeV/fm^. Further- 
more it is difficult to adjust the latent heat independently 
of Tc in such models of the phase transition. 

We have taken a more pragmatic approach and have 
constructed a thermodynamically consistent EOS with a 
variable latent heat in the e and Hb plane. First, note 
the following two derivatives which apply along the path 
where ub/s is constant 



ub/s 



TLB I S 



p + e 



(10) 
(11) 



The first of these is simply the definition of the speed 
of sound. The second relation is surprising: it does not 
contain the chemical potential /is explicitly. (It follows 

by noting that (^) , ^ (¥) + ( :f^) — (tt) , 

o \ del nB I s \de/7iB ydriB ) ^ s \de I tlb I s 

and solving for ^ /s' ^-^ using thermodynamic iden- 
tities). Given the speed of sound everywhere and the en- 
tropy on a single arc in the e, plane, these derivatives 
may be integrated to determine the entropy, s(e,nB). 
From the entropy, all other thermodynamic functions 
(e.g., T and \ib^ may be determined. Below, only the 
speed of sound is specified. 

For smooth flows, entropy and baryon number are sep- 
arately conserved. If at some initial time ub / s is constant 
everywhere in space, the two conservation laws i mply 
that ubI s is constant everywhere in space and time [ [40[ . 
For the initial conditions specified in Sect. II B, ub/s — 



Cub/Cs is constant in space and remains constant as the 
system evolves. Therefore, the pressure is needed only 
along the path ub/s = Cn-^jCs- This may be directly 
verified by fully differentiating — and noting 

that the derivatives of the pressure only appear as the 

speedofsound,(^)^^^^.c-(if)^^4-f-(^)^. 

Strictly speaking, transverse shock waves develop near 
the phase transition and invalidate the assumption of en- 
tropy conservation. However, numerical and analytical 
evidence has shown that entropy production in hydro- 
dynamic simulations of nucleus-nucleus collisions is at 
most a few percent Below, entropy production is 

ignored and the pressure is specified along the trajectory 

The EOS consists of three pieces: a hadronic phase, 
a mixed phase, and a QGP phase. In strong interac- 
tions, Baryon number (B), Strangeness (S), and Isospin 
(I) are conserved and therefore the EOS depends on T 
and /iB,Ais,and /i/. In the hadronic phase, the thermo- 
dynamic quantities -the pressure (p), the energy den- 
sity (e), the entropy density (s), and number densities 
(jiQ where Q=B,S,I)- are taken as ideal gas mixtures of 
the lowest SU(3) multiplets of mesons and baryons. The 



mix includes the pseudo-scalar meson octet (tt, ry, K) and 
singlet (?]'), the vector meson octet {p,K*,uj) and sin- 
glet (0), the i baryon and anti-baryon octets and the 

I baryon and the anti-baryon decuplets. 
p, e, s and uq are given by 



Specifically, 



n-Q ^ (T, p.BBi + fisSi + fiih) (12) 

i 

P = ^Pid A^B-Bi + fisSi + /i//i) (13) 

i 

e = '^e^^{T,fiBBt + p,sS.t + pili) (14) 

i 

s = ^s'^^{T,^bBi + p.sSi+ fiili) . (15) 

i 

Here the sum is over the hadrons species, B^^Si^Ii are 
the quantum numbers of the i-th hadron, Ui is + for 
bosons but — for fermions, and for example, p^(r, /i) is 
the pressure of a simple ideal Bose gas. A fast numerical 
method for evaluating the thermodynamic quantities of 
simple Bose/Fermi gases has been constructed ||5j]. For a 
given T, (/is, /i^, /i/) are determined by the requirements 
that total strangeness {ns) and isospin(n/) be zero and 
that tib/s = CubI^^s- The thermodynamic quantities 
are then taken as functions of e along the adiabatic path 
where ub/S = C„g/Cs. This hadronic EOS is taken up 
to a temperature of Tc = 165 MeV or an energy density 
en ~ 0.45 GeV/fm'^ (see Fig.||). The squared speed of 
sound is approximately 1/5 in this hadronic gas. 

Above the hadronic phase, only the speed of sound 
squared, c^, is specified. For the mixed phase the speed 
of sound was made approximately zero, = 0.02 c. The 
width of the mixed phase (see Fig. ^) is the Latent Heat 
(LH), LH= eg — en- LH is taken as a parameter and 
is adjusted to form phase diagrams LH8, LH16,... with 
latent heats, 0.80 GeV/fm^ 1.6 GeV/fm^.. . Above the 
mixed phase, e > eg, the degrees of freedom are taken 
as massless and the speed of sound is accordingly, Cs — 
■\/l/3. We also consider two limiting cases: a Resonance 
Gas (RG) EOS and LUoo. For a RG EOS, the speed of 
sound is constant above en- For LHoo, the mixed phase 
continues forever [eg = oo) and there is no ideal plasma 
phase. 

With the speed of sound specified in all phases, Eq.p^ 
and ^ are integrated to find the pressure and entropy 
along the adiabatic path specified by the initial condi- 
tions, Ub/s — Cnsl^s- The full phase diagram for SPS 
initial conditions is shown in Fig.|^. In Sect. IV and 
Sect. ^ the subset of the EOS consistent with the avail- 
able radial and elliptic flow data is found. 



D. RQMD and the Cooper Frye Formula 

Given the initial conditions and the EOS, the equations 
of motion are integrated in time. As the system expands 
and cools, the mean free path becomes much less than 
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persurface and f^{E,iJ,) is the Bose/Fermi distribution 




function 



a. I 



e(GeV/fm ) 



FIG. 2: The pressure (p) versus energy density (e) for dif- 
ferent EOS. EOS LH8, LH16 and LHoo become increasingly 
soft and have latent heats 0.8GeV/fm"^, 1.6GeV/fm"^, and 
oo. The EOS are shown along the adiabatic path for SPS 
initial conditions, s/ub ~ 42. For RHIC initial conditions, 
s/hb ~ 150, the changes are small. 



the nuclear radius, and the system breaks up into free- 
streaming particles. Typically, hydro practitioners fl^ , 
[iTt model the breakup or freezeout stage by finding a 
hypersurface in space and time where the temperature 
equals some freezeout temperature, T = Tj. 

This simple picture in which all particles are emitted 
from a single space-time surface is however unrealistic. 
Different particles have different hadronic cross sections 
and suffer their last interaction at different times. They 
are emitted over a space-time region rather than on a 
sharp surface. Further, particles in the center rescatter 
for a longer time than particles in the periphery. To 
model this physics, the spectrum of particles exiting the 
space-time surface is taken as the input to the hadronic 
cascade, RQMD 
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Subsequently, the particles rescat- 
ter. Below, the space-time surface is referred to as the 
switching surface rather than the freezeout surface. 
The attending problems with this approach are described 
below after the details of the input distribution to RQMD 
are described. 

For the family of EOS discussed above, for T < Tc, the 
fluid is made up of a collection of ideal gases of fermions 
and bosons. In order to conserve energy and momentum 
across the surface, the spectrum of (species) i crossing the 
surface is given by 139, M , 



(Pp 



fid (P ■ f^sBi + fisSi + fiili) 

X dS„ (16) 



The particle index i runs over all 

the species in the EOS - no more, no less. If the quasi- 
particles are interacting and the EOS is non-ideal, due 
to viscosity, mean fields, particle lifetimes, etc., then the 
should be modified accordingly. 
The differential elements of the hypersurface d'Ea can 
be separated into time-like (dSadS" > 0) and space-like 
(dEo-dS" < 0) surface elements. For time-like surfaces, 
the integrand in Eq.|l^ is positive and there is a frame 
(the rest frame of the surface), where dJ^a = (dV, 0,0,0). 
The spectrum of Eq. ^ is simply a thermal spectrum 
boosted by the flow velocity in the frame of the surface. 
(In practice, the surface velocity is small for most time- 
like surfaces). The yield of (species) ^ leaving a surface el- 
ement, dT,a, is simply n^^(T, fj.BBi + fj.sSi + fJ.iIi)U°'dT,a, 
as may be found by integrating the left and right sides 

,3 

with and going to the rest frame of the matter. 

For space-like surfaces, the integrand in Eq.|l^ is both 
positive and negative depending on the momentum of the 
particle. When the integrand is positive, the particle is 
leaving the surface and when it is negative the particle is 
entering the surface. We reject particles entering the hy- 
drodynamic surface. The distribution of particles exiting 
a space-like surface is 



d^N- 



d^p 



X dSc,e(p"dS]a) . 



(17) 



where dYl^ is a differential element of the freezeout hy- 



It is this distribution that we take as the input distri- 
bution for RQMD. For a discussion of the problem of 
space-like surfaces see The number of particles 

leaving a differential surface element is given by a more 
complicated formula which is again found by integrating 
both sides of the equation with ^ ||5^ . For a stationary 
surface in the rest frame, the formula has the simple in- 
terpretation as the number of particles evaporating from 
the surface per area per unit time. A consequence of the 
theta-function is that energy, momentum, and particle 
number are not exactly conserved across the transition 
surface. However, the error is only « 2.0%(« 5.0%) for 
central (peripheral b=8.0fm) AuAu collisions at RHIC. 



III. THE SPACE-TIME EVOLUTION 
A. The Hydrodynamic Solution 

This section reviews the hydrodynamic evolution for 
different EOS used in this work. The evolution at the 
SPS and RHIC is summarized in Fig. ||. The switch- 
ing isotherm, en (shown in the middle), is particularly 
important since in the hydro-|-cascade approach, the par- 
ticles are injected into RQMD with the velocity distribu- 
tion of this isotherm. 

For EOS with a phase transition (LH8) there are three 
phases and three corresponding stages in the acceleration 
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FIG. 3: The left and right sides show the hydrodynamic 
solution at the SPS and RHIC for different EOS. The thin 
lines show contours of constant transverse fluid rapidity {vt ~ 
tanh(j/T)) with values 0.1,0.2,.. .,0.7 . The thick lines show 
contours of constant energy density. ei2o denotes the energy 
density where T = 120 MeV. ch and eg denote the energy 
density (for LH8) where the matter shifts from hadronic to 
mixed and mixed to a QGP, respectively. The shift to RQMD 
is made at eh- {vt) denotes the mean transverse rapidity 
weighted with the total entropy flowing through the energy 
density contours. Walking along these contours, the line is 
broken into segments by dashed and then solid lines. 20% 
of the total entropy passing through the entire arc passes 
through each segment. 



history, (i) an explosive QGP phase (e > eg), in which 
the matter accelerates rapidly, (ii) a soft mixed phase 
(en < e < eg), in which the matter free streams with 
constant velocity and (iii) a hadronic phase (e < en), in 
which the hadronic pressure produces additional acceler- 
ation. 

The QGP phase dictates the duration and transverse 
size of the collision. At RHIC, the QGP pressure drives 
the matter outward, rapidly increasing the radius, which 
in turn shortens the overall lifetime. Therefore, approx- 



imately doubling the total multiplicity from the SPS 
to RHIC increases the total lifetime only slightly, from 
lOfm/c to 11 fm/c. All the additional multiplicity is ab- 
sorbed by a slightly larger transverse radius. Similarly, 
for a RG EOS the acceleration is robust and continu- 
ous and increasing the total multiplicity only slightly in- 
creases the radius and lifetime. In bulk, the radii and 
lifetimes of RG are similar to LH8. 

By contrast, for LHoo, the stiff QGP phase is absent 
and the mixed phase is dominant at high energy densi- 
ties. The strong transverse acceleration associated with 
LH8, is replaced with a slow evaporative process. The 
radius of the system slowly shrinks as a function of time. 
Unlike LH8, increasing the total multiplicity increases 
the lifetime rather than the radius. Between the SPS 
and RHIC the lifetime increases from from 14 fm/c to 
21 fm/c. Summarizing, the QGP drives a transverse ex- 
pansion; the transverse expansion increases the radius 
and shortens the overall lifetime compared to an EOS 
without the QGP push. 

To quantify the input velocity distributions into 
RQMD, Fig. ^ plots the transverse fluid rapidity ver- 
sus the radius along the switching isotherm en at the 
SPS and RHIC. For LH8 and RG, the transverse rapid- 
ity shows a linear rise with radius. A linear flow profile is 
often used in phenomcnological fits to the particle spectra 
; this calculation validates this approach. 

For an EOS with a phase transition to the QGP (LH8), 
the acceleration is initially large but subsequently stalls 
in the mixed phase. By contrast, for an EOS without 
the phase transition (RG), the acceleration is robust and 
continuous. Therefore, although the initial transverse ac- 
celeration is smaller for a RG than for LH8, the RG veloc- 
ity at the end of the SPS mixed phase is 50% larger than 
for LH8. At RHIC, where for LH8 the QGP phase lives 
longer, the RG velocity is only « 15% larger. Although 
mean lifetimes and radii of the RG EOS are similar to 
LH8, the change in the velocity distributions from the 
SPS to RHIC are markedly different. Comparing LH8 
to LHoo at RHIC, the LH8 flow velocity is approximately 
twice as large as the LHoo flow velocity. 

Nevertheless, it should be noted that if freezeout is 
taken as Tf « 120 McV, then the differences between the 
flow velocities of the EOS is smeared out by the hadron 
phase, as can be seen by examining the mean flow veloc- 
ities on the ei2o curves in Fig.|^. Indeed, the hadronic 
phase of LHoo (which lives longer, since it is born with 
no transverse velocity) can partially compensate for the 
weak initial acceleration. Along the Tf « 120 MeV 
isotherm, the flow velocities of LHoo, RG and LH8 are 
roughly comparable. 

To characterize the flow in non-central collisions for 
the EOS used in this work, we follow Kolb et al. jl^ and 
calculate a quantity derived from the stress tensor for 



(18) 



where the ()s denotes an average over the transverse 
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FIG. 4: Walking along the eh contours in Fig. ^ (where the 
switch to RQMD is made), the transverse rapidity is traced 
as a function of radius at (a) the SPS and (b) RHIC for the 
three EOS. See Fig. ^ for an explanation of the dashed and 
solid segments. 



plane weighted with the entropy per area per unit spatial 

the final mo- 



rapidity, S7 rdx dy. Cp is related to 



mentum anisotropy of the particle distribution [g9| which 
ultimately is related to V2{pt) @|- 

The general trends seen in ep follow from the discussion 
above on the hydrodynamic response of each EOS. LH8 
shows a strong early response followed by a stall and 
subsequent flattening as the matter distribution becomes 
almost spherical. At RHIC the strong early response 
lives substantially longer as the matter spends a larger 
fraction of its total lifetime in the QGP phase. For LHoo 
the matter is initially stalled but rapidly accelerates as it 
slowly enters the hadronic phase (see Fig. (e) and (f)). 
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FIG. 5: Following Kolb et al. we show the anisotropy of 
the stress tensor, ep (a measure of elliptic flow) for the EOS 
used in this work. The solid (dashed) curves are for RHIC 
(SPS). The solid (open) symbols indicate when the center of 
the fluid passes through a temperature of 160 MeV (120 MeV). 
The solid symbols are therefore representative of the switching 
temperature to RQMD or en in Fig. ^ 



B. Qualitative Predictions of the Hydrodynamic 
Response 

We can now make some qualitative predictions from 
the hydrodynamic solution. Assume momentarily that 
the final hadron momentum distributions reflect the 
boundary between the mixed and hadronic phases or ch 
in Fig. ||. Then with the curves presented in the last sec- 
tion, LH8 predicts two qualitative changes. First with 
Fig. |, between the SPS and RHIC, weighted ellip- 
tic flow should increase by almost a factor of two as the 
QGP replaces the mixed phase and dominates the early 
evolution. Second with Fig. ^, the total transverse mo- 
mentum should increase by 30% as the QGP drives ad- 
ditional transverse motion. The flow differentiates LH8 
from a RG EOS and from LHoo. Since RG EOS acceler- 
ates continuously and does not stall in the mixed phase, 
the transverse momentum is larger at both the SPS and 
RHIC. In addition, for SPS collision energies, the elliptic 
flow (cp) for a RG EOS is almost a factor of two larger 
than for LH8. For LHcxd, the transverse momentum is 
very low until the very end. 

To make these qualitative predictions quantitative and 
to compare the hydrodynamic solution to experimental 
data, it is essential to model hadronic freezeout. Between 
the time when Tf = 160 MeV (the solid symbols) and 
Tf — 120 MeV (the open symbols) elliptic flow changes 
dramatically for each EOS. This is especially true for 
LH8 at the SPS, where the mixed phase abruptly stalls 
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FIG. 6; The fraction of pions emitted into RQMD per unit 
of time relative to the last instant (rmai) that the matter is 
evolved by hydrodynamics. 



the development of elliptic flow but then the hadronic 
phase rapidly completes the development. The differ- 
ences in the early acceleration tend to get washed out 
by the hadronic stage. Indeed, even LHcx) has a rea- 
sonable radial and elliptic flow by Tf = 120 MeV. The 
extent to which signatures of the early QGP accelera- 
tion remain in the final spectra depends on whether the 
freezeout temperature should be taken as Tf = 120 MeV 
or Tf — 160 MeV. The breakup of a heavy ion collision 
can only be addressed with hadronic cross sections and 
expansion rates. 



C. RQMD — Input and Response 

Relativistic Quantum Molecular Dynamics (RQMD) 
p6[ is a hadronic transport computer code which incor- 
porates many known hadronic cross sections. RQMD 
has been used extensively to model the heavy ion dy- 
namics |o[ ^ |2|, Briefly, when two particles 
come within d < -^/ct/tt, they elastically scatter or form a 
resonance. Resonance formation and decay dominate the 
evolution. The principal reactions are tttt — > p, tt iV ^ A, 
TT K ^ K*, and ttA — > E*. Only binary collisions are 
considered in this hadron cascade. Before discussing the 
response of RQMD, we first consider the input. 

The time distributions are found by projecting the en- 
tropy distribution on the switching isotherm ch in Fig. ^ 
(or T = 160 MeV) onto the r axis. In Fig. || the frac- 
tion of pions (or entropy) injected into RQMD per unit 
r is plotted as a function r/Tmax for each EOS. For LH8, 
very few particles are evaporated from space-like surfaces 



at early times, and at r/Tmax ~ 0.8 (t « 9fm/c) parti- 
cles are emitted in bulk from the time component of the 
switching surface. For LHcxd, particles are continuously 
evaporated from the transition surface and the radius 
slowly decreases. Therefore, the time distribution is rel- 
atively uniform. Finally for a RG EOS, the freezeout 
surface is not box-like and particles are also emitted into 
RQMD slowly and continuously. 

Now consider the dynamic response of the hadronic 
cascade. For LH8, the hydrodynamic input into RQMD 
can be characterized as a simple thermal model with a 
linearly rising flow profile with a uniform radial distribu- 
tion except at the edge of the distribution where there is 
a small maximum. Once this input distribution is taken, 
the hadrons re-scatter within RQMD and different parti- 
cles decouple from the cascade at different times. Fig. |^ 
(a) plots the mean emission time < r > (the time of 
last interaction) versus the mass of the particle species. 
Also shown is < r >, when all collisions in RQMD are 
switched off and only resonance decays are allowed. The 
mean number of collisions experienced by a particle is 
shown in Fig. |^ (b). The mesons scatter approximately 
once after their principal resonances (p, i4r*,etc.) decay 
and decouple around r « 14fm/c. In contrast, due to 
strong meson-baryon resonances A, S*, ... , nucleons and 
hyperons (A and S) scatter approximately twice and de- 
couple around r « 18fm/c. The <j) and are emitted 
directly from the phase boundary since they have small 
hadronic cross sections. 

The duration of the hadronic stage dictates the spatial 
extent of the final source. In Fig |^ (c), the mean radius 
is shown as a function of particle mass with and without 
re-scattering in the hadronic cascade. For comparison, 
we apply the simple formulas: We assume all particles 
are emitted from the switching surface at a mean radius 
Ro and a mean time Tq, with a constant radial velocity vt 
(see Fig. Q (a) and (c)). Since 4> is emitted directly from 
the switching surface, Ro ={R,p) and Tq = (r^). With 
the formula distance = velocity x time, we have 



(R)^ = Ro + Vt{{Tx) - To), 



(19) 



where {R)x {t^) is the freezeout radius (time) of particle 
X and Vt is the freezeout drift velocity. This velocity 
incorporates a thermal drift velocity and the flow velocity 
of the source. Given a constant velocity as a function of 
mass Vt ~ .42 c, a very simple fit to the freezeout radii 
is obtained, as shown in Fig. 1(c). Thus, hadr onic cross 
sections dictate the final freezeout radii of the source. 

Hadronic cross-sections also dictate the spatial geome- 
try in non-central collisions. In non-central collisions the 
ellipticity of the source at freezeout is quantified by the 
spatial anisotropy, 



S2 



y 



(20) 



Here, the averages are taken over points of last interac- 
tion in the cascade. S2 is negative for the initial almond- 
shaped distribution but positive for a cucumber-shaped 
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FIG. 7: The mean (a) emission time, (b) number of collisions, 
(c) emission radius, and (d) spatial anisotropy S2 (see Eq. pc| ) 
as a function of particle mass with and without the RQMD 
hadronic after-burner. The averages are taken over the points 
of last interaction for AuAu collisions at b=6fm at RHIC. In 
(a) and (c). To, and Ro label the mean emission time and 
radius of the meson, vt denotes the freezeout drift velocity 
of Eq. 



distribution. Fig. |^(d) shows S2 without re-scattering but 
with resonance decays (Hydro Only) and with hadronic 
re-scattering (Hydro-|-RQMD). The initial elliptic flow 
(^2) changes the overall geometry (52) by the end of the 
RQMD stage. At the end of the hydrodynamic stage 
S2 is negative, indicating that the source retains at least 
some of its initial almond distribution. S2 becomes posi- 
tive as the system evolves and the momentum asymmetry 
changes the source geometry. For nucleons, S2 is almost 
for modest impact parameters; this may have ob- 
servable consequences p3] . 



D. Impact Parameter Dependence of the Space 
Time Evolution 

In the previous section, we discussed how hadronic 
cross sections control the lifetime and geometry of the 
final hadronic distributions. Now the impact parame- 
ter is varied and the freezeout distributions are modified. 
In non-central collisions, the number of charged parti- 
cles scales as the number of participants; therefore the 
lifetime of the hadronic stage should also scale as the 
number of participants. However, the hadronic lifetime 
is also a function of the cross section, the radius, and 
the expansion rate (df^U^). These depend respectively 
on the particle species, the r.m.s. radius of the initial 
geometry, and the EOS. In Fig. ^, the different contribu- 
tions to the total lifetime are studied. We plot the mean 
emission time < t > , divided by size R^J^^'"^^ , as a func- 
tion of the number of participants for different particles 
and EOS. To set the absolute scale, the "free" axes show 
< r > directly at two impact parameters. 

Consider first the LH8 curves (a): The total lifetime 
for all particle species falls by approximately 30% from 
central (b=Ofm) to peripheral (b=8fm) collisions. The 
order of particle emission remains as the impact param- 
eter is varied: First rare species (</), Q) are emitted, then 
mesons (tt, if) and finally baryons (A^, A). For the 0, 
which is representative of the hydrodynamic stage, the 
curves in Fig. ^(a) are flat at the 15% level, indicating 
that the total lifetime scales roughly with the size of the 
overlap region. For pions, the total lifetime also scales 
with Rrms- For protons, indicative of baryon emission, 
the total lifetime does not quite scale as Rrms but rather 
depends on the absolute number of charged particles in 
addition to the geometry. This is natural since the freeze- 
out of protons is controlled by the formation of A reso- 
nances. 

For LHc», (r) /Rrms does change more rapidly than for 
LH8 This is especially true for nucleons. However, for cf) 
and TT the difference in the Np dependence of (t) /Rrms 
is small and to a reasonable approximation, the lifetimes 
of and TT scale with R'^/^f"'' for all EOS. Changing the 
EOS simply moves the various curves up and down on 
Fig. H (a) and (b). A RG EOS was also studied (not 
shown) and the lifetime and 
similar to LH8. 



Np dependence were quite 
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FIG. 8: The mean emission time (r) relative to the Glauber 
r.m.s. radius (see Eq. ^) as a function of the number of par- 
ticipating nucleons for (a) LH8 and (b) LHoo. The free axes 
show the absolute lifetime at impact parameters b=Ofm and 
b=8 fm. 



Eq. |l^ relates the freezeout radii and geometry of the 
different particles to a freezeout time and a single freeze- 
out drift velocity (vt) and a single freezeout radius (i?^). 
This simple formula was found to be applicable to all im- 
pact parameters with approximately the same drift veloc- 
ity as in central collisions. The lifetimes and radii (i?^) 
all scale with the root mean square radius. 

Given the rather simple scaling of lifetimes and radii as 
a function of impact parameter, it is natural to consider 
the density of pions at freezeout as first done in 
Since pion number is approximately conserved during the 
cascading process, we have s oc rt^ and the freezeout 



FIG. 9: The density of pions at freezeout (see text) versus 
the number of participants for different EOS. Here < Tn > 
and < > denote the mean pion emission time and radius. 



entropy density is 



dN^/dy 



n < T >< R >2 



(21) 



This quantity is shown as a function of the number of 
participants for RHIC collisions in Fig. ^. The freeze- 
out entropy density is roughly constant as a function of 
impact parameter. In addition, the freezeout density is 
independent of EOS in spite of differences in transverse 
velocity gradients. It has been argued that central PbPb 
collisions cool to a lower temperature than peripheral 
collisions since transverse and longitudinal velocity gra- 
dients are larger in peripheral collisions jBSf . However 
at least in the model, freezeout is not driven by the ex- 
pansion rate; rather, the freezeout condition reflects a 
density where the mean free path becomes comparable 
to the radius of the nucleus. 

Next we hold the collision geometry fixed and examine 
the changes in lifetime, radius and freezeout density as 
the initial entropy density (i.e., the collision multiplicity) 
is increased. Fig. 
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(a) and (b) shows the lifetime and 
emission radius for a PbPb collision at b=6 fm as the 
multiplicity per participant is increased from the SPS to 
the RHIC domain. For LH8, the radius increases but the 
emission time does not, while for LHoo the situation is 
exactly reversed. 

This behavior is readily understood as entropy con- 
servation in the transverse plane. Entropy conservation 
(see Eq. ^ relates the entropy density to the Bjorken 
time (r), the total conserved entropy per unit rapidity 
(dStot/dy), and the effective area (Aeff) of the source 
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FIG. 10: The mean pion (a) emission time, (b) emission 
radius, and (c) freezeout density as a function of the total 
muhiplicity in a PbPb coUision at b=6fm. The SPS and 
RHIC arrows denote where the charged particle per partici- 
pant per rapidity matches SPS and RHIC initial conditions 
respectively. 



with the schematic relation, 

dStot/dy 



< S >r 



(22) 



As seen with Fig.^(c), the entropy density at freezeout 



(s/) is roughly constant as a function of beam energy. 
The freezeout time Tf may therefore be related to freeze- 
out entropy density s / with 



dStot/dy 



(23) 



For LH8, the strong transverse acceleration rapidly in- 
creases the area and lowers the entropy density s to s/. 
Consequently, as the multiplicity is doubled, the total 
lifetime increases by only 20%. For LHoo, the radius does 
not increase but the lifetime increases significantly. Thus 
the transverse expansion, together with entropy conser- 
vation, ultimately determine the total lifetime. 



IV. RADIAL FLOW FROM THE SPS TO RHIC 
A. The SPS 

In a traditional hydrodynamic calculation ^ |5j , 
the pion and nucleon yields fix the total entropy and 
baryon number in the initial conditions. The freezeout 
temperature is adjusted to fit the pion and proton pT 
spectra. In the hydro-|-cascade approach advocated here, 
the freezeout temperature is not a parameter since parti- 
cles decouple from the cascade when their collision rates 
become small. Therefore, the pion and nucleon yields set 
the total entropy and baryon number as before, but the 
slope parameters provide significant information about 
the EOS. In particular, the latent heat of the phase tran- 
sition which best matches the pion and nucleon pT spec- 
trum is LH w 0.8 GeV/fm^ 

In the previous section, we discussed a family of EOS, 
each with a different latent heat. Now we show in Fig. |l|, 
the calculated net proton spectrum for the resonance gas 
EOS and for EOS LH4-LH16. The experimental and the- 
oretical spectra are absolutely normalized. The two num- 
bers parameterizing the initial conditions Cs and C„g 
(see Table |) are adjusted to match the height of these 
spectra. The model curves have been multiplied by a fac- 
tor of 0.93 to account for the fact that the data is 5% cen- 
tral. Once the height of the spectrum is tuned, the shape 
of the spectrum is determined by the course of the hydro- 
dynamic evolution, or more generally, by pressure gradi- 
ents and the duration of the collision. Therefore, it is sig- 
nificant that hydrodynamics generates a flow vt ~ 0.5, 
which is needed to explain the spectra. This flow velocity 
was extracted from a variety of thermal analyses [l^ . 

For EOS with large latent heats (e.g. LH16), the 
Pt spectrum is too soft. This is because the hydrody- 
namic system spends a long time in the mixed phase in 
which pressure gradients do not generate collective mo- 
tion. Bag-model equations of state, employed in many 
hydrodynamic calculations ^2] , typically have a la- 
tent heat from 1 — 1.5GeV/fm which makes the EOS 
rather soft. This large latent heat is usuall y c ompensated 
by adjusting the freezeout temperature |1^. An EOS 
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FIG. 11: A comparison to NA49 (a) negative hadron 
and (b) net-proton spectra for different EOS. 



with a modest first order phase transition (e.g. LH4 and 
LH8) generally reproduces the shape of the pr spectra in 
FigO. Unfortunately, a RG EOS can also reproduce the 
shape of the px spectrum and additional experimental 
information is needed to separate EOS. 

The slope systematics of A, S and n provide the nec- 
essary information. The presence of a phase transi- 
tion stalls the acceleration |2^, therefore informa- 
tion about the velocity at the end of the mixed phase 
can separate a RG EOS from LH8. At the SPS, the 
spectra of the different particle species are all reasonably 
exponential and a slope parameter, T^iope is extracted. 



1 dN 
Mt dMr 



— Cexp(- 



slope 



-), 



(24) 



■ Pr- 



ill this parameterization, 
is directly related to the 
The model spectra are fit 



where Mt — \/ vn? 
the slope parameter Tslope^ 

< Mt >= 2Tslope + T,,",-Hm- 

with Eq. ^ over the range corresponding to the WA97 
[ ]68[ experimental acceptance (Mt — m = 0.0 — 0.9 GeV); 
the slope parameters are shown versus particle mass with 
and without RQMD in Fig. |l|(a). First look at the "Hy- 
dro Only" curves. The slopes increase approximately lin- 
early with mass as is expected in a thermal expanding 
source model ||l^, |5^, ^ . The non-monotonous increas- 
ing mass in the "Hydro Only" curves is due to resonance 
decays and the baryon content of the particles. Once 
RQMD is included, the slopes are modified by hadronic 
rescattering leading to mass dependence characteristic of 

Note the following features. 



differential freezeout |34 



First, in the model the gives a good measure of the fiow 
velocity at the end of the mixed phase. Second, note the 
« 40% increase in the nucleon and A slope parameters 
and the small decrease in the pion slope parameter, due 
to cooling. As the hadron gas expands, the pions excite 
A and S* resonances and drive additional transverse mo- 
tion in the nucleon and hyperon sectors. However, the 
pions increase the nucleon (Mt) only at the expense of 
their own kinetic energy. In a traditional hydrodynamic 
approach, the hydrodynamic evolution would be contin- 
ued to match the slope of the nucleon spectrum. Judging 
from Fig. ^(a), this is misguided, as a nucleon receives 
much of its momentum after pions have decoupled. As 
shown below, the nucleon receives about 20% of its trans- 
verse kinetic energy from the pion "wind" , irrespective 
of the colliding energy. To incorporate the rich freezeout 
dynamics of a cascade, different freezeout temperatures 
and velocities should be taken for different particles . 

A comparison to the available data on slope systemat- 
ics is given in Fig.|l^: (b) shows the slopes for different 
types of EOS while (c) shows the sensitivity to the la- 
tent heat. Although RG and LH4 are capable of repro- 
ducing the pion and nucleon spectra, they significantly 
over-predict the slope parameters of A, S and 51^. This 
is because LH4 and RG already have developed a sub- 
stantial fiow velocity at the end of the mixed phase. The 
slope parameter of the is a sensitive measure of the 
flow velocity at the end of the mixed phase since the 
flow velocity is amplified by the mass in the approximate 
formula Tgiope = Ttt + m < v >^. LHoo, by contrast, 
under-predicts the slope parameters of A and S, indicat- 
ing that (vt) is too small at the end of the mixed phase. 
With LH8, the acceleration is modest - but significant - 
and the slope systematics are generally reproduced. In 
the model, only an EOS which has both a stiff and soft 
piece is capable of reproducing general trends seen in the 
particle spectra. 
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FIG. 13: The nucleon Mr spectrum at the SPS and RHIC 
with and without the RQMD after-burner. 



B. Qualitative Changes at RHIC 

It was argued above that LH8 provides the best de- 
scription of the radial flow at SPS: Now the same EOS is 
used to make predictions for RHIC. At RHIC, the initial 
energy density is well above the phase transition, and the 
large early pressure is expected to drive collective motion. 
In Fig. |l^, the nucleon Mt spectrum for the SPS and 
RHIC are shown with and without the hadronic rescat- 
tering in RQMD. Two features are immediately observed: 
1. The (Mt) increases as the collision energy is increased 
from the SPS to RHIC H, ||. 2. The spectra with- 
out hadronic rescattering are reasonably well described 
by a single exponential (i.e., they look linear on the log 
plot shown). Once rescattering is included the spectra 
are curved; this curvature grows from the SPS to RHIC. 
Describing the spectra with a single slope parameter, al- 
though useful in summarizing a large variety of data, is 
only approximate. 



C. The (Mt) from the SPS to RHIC and Beyond 



FIG. 12: A compilation of slope parameters from the SPS 
experiments and a comparison to the model for diflferent 
EOS. (a) shows the model slope parameters with and with- 
out hadronic rescattering (b) shows the slope parameters for 
three qualitatively different EOS: RG, LH8, and LHoo. (c) 
shows the sensitivity to the latent heat in the LH family of 
EOS. NA49 data points are from ||] while WA97 data points 
are from j68| . The model spectra are fit with Eq. ^ over the 
range corresponding to the WA97 experimental accep- 
tance (Mt - m = 0.0 - 0.9 GeV). 



To summarize the bulk energy transport in the model 
we show in Fig.|l^: (a) the (Mt) — M for nucleons as a 
function of collision energy (dNch/dy) for different EOS, 
and (b) the (Mt) — M for different particle species as a 
function of mass for different collision energies (particle 
multiplicities). 

Fig. |l^(a) demonstrates that the additional entropy 
gets converted into additional transverse motion for each 



15 



0.6 



> 0.5 
(5 



A 
I- 



0.4 



0.3 



0.2 



* RG EoS 


,-■* 


^LH8 EoS 




• LHoo 


EoS ^ 








,.■* 


_ *■■' 


/ ■■* 


• 


* 






■ .-■■;'-•'* 






~ *"' 








.«■■" 












40(^ 


600 


SPS 


RHIC 


dN,,/dy 



> 

0) 

O 



A 
I- 

V 



0.7 



0.6 



0.5 



0.4 



0.3 



0.2 



dN/dy = 220. 
dN/dy = 345 
dN/dy = 466 
dN/dy = 699 




AZ S 
■f + + 



n 



0.2 0.4 0.6 0.8 1 



1.2 1.4 1.6 1.8 
M(GeV) 



FIG. 14: For PbPb collisions at b=6fm, (a) shows 
{Mt)—M for nucleons versus the charged particle multiplicity 
{dNch/dy). For each EOS the upper curve is with hadronic 
rescattering (Hydro+RQMD) and the lower curve is without 
hadronic rescattering (Hydro Only), (b) shows (A/t) — M for 
different particle species as the total multiplicity is increased 
from the SPS to RHIC and beyond. 



EOS. For each EOS, the hadronic contribution to the 
mean (Mr) remains constant and is approximately 20% 
for LH8 and RG, but is approximately 30% for LHoo. 
LHoo is a special case, and we may say that the RQMD 
contribution is approximately 20% and is independent of 
the underlying EOS. 

Fig. |lj(b) demonstrates how the increase in the mean 
Mt influences the spectra of different particles by plot- 
ting {Mt) versus mass iQ. At the SPS the flow veloc- 
ity at the end of the mixed phase is relatively small - 
vt ~ 0.4 . The slopes before the RQMD phase show a 



linear rise characteristic of hydro, Tsiope = Tth + m (w^). 
When the flow velocity is small, hadronic rescattering 
changes the linear mass dependence significantly, giving 
the characteristic shape observed at the SPS. As the flow 
velocity increases from the SPS to RHIC and beyond, 
the linear rise with mass becomes increasingly steep and 
hadronic rescattering, while still contributing to 20% of 
the (Mt) for nucleons, does not change the overall mass 
dependence. The qualitative shape of the mass depen- 
dence of (Mt) therefore gives a good measure of the flow 
velocity at the end of the mixed phase. Since this flow 
velocity is different for different EOS, the mass depen- 
dence of the (Mt) can therefore separate the different 
EOS studied. 



D. The Flow Profile from the SPS to RHIC 

The curvature in the Mt spectrum is a signature of a 
radially flowing thermal source. The general features can 
be understood from a simple thermal model. For a cylin- 
drically symmetric shell, which expands longitudinally in 
a boost invariant fashion and which freezes out in an in- 
stant with constant temperature T, and a radial velocity 
Vt = tanhp, the Mt spectrum is given by ^ 



1 



WtIm-t ^ ^^^Th[—^]K, 



For pt ^ m we have 



Mt cosh p 



Tt 



th 



- slope 



th 



1 



(25) 



(26) 



Generally, increasing the velocity increases the curva- 
ture of the final spectrum for heavy particles. Increasing 
the mass also increases the curvature of the flow profile. 
The shape of the spectrum, together with the mass de- 
pendence of the observed particle, may provide a good 
experimental measure of the velocity of the source at 
hadronization. 

Hadronic rescattering changes the curvature seen in 
the thermal spectra discussed above. First, different 
particle types have different hadronic cross sections and 
therefore freezeout at different times and with different 
velocities. The curvatures in the final spectra measure 
these different freezeout velocities. Second, the cascade 
generates additional transverse flow predominantly in the 
low Mt region of the spectra. To quantify these ef- 
fects, we divide the Mt spectra into a low Mt region, 
< Mt < 0.6 GeV and a high Mt region, 0.6 GeV < 
Mt < 1.6 GeV. We then fit an exponential in both do- 
mains. Thus, there is a low Mt slope and high Mt slope. 
We have checked that this parameterization gives a good 
description of the shape for all the spectra considered. 
Fig. |5| shows the low and high Mt slopes as a function 
of the particle mass, with and without the RQMD after- 
burner. The curves illustrate the collective acceleration 
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FIG. 15: Each particle spectrum in central AuAu collisions, 
MtMIt ' with an exponential slope parameter from < 
Mt <^0.6GeV and from 0.6 < Mr < 1.6 GeV. Thus there 
is a low and a high A/t slope representative of the slopes 
measured by the STAR Q and PHENIX collaborations @. 
(a) shows the mass dependence of the low and high Mt slopes 
with and without the RQMD after-burner for LH8. (b) shows 
the mass dependence of the low and high Mt slopes for the 
different EOS with the RQMD after burner. 



which occurs during hadronic rescattering, and iUustrate 
an interplay between freezeout and hydrodynamic behav- 
ior. First look at the "Hydro Only" curves: the curva- 
ture (i.e. the difference between the low Mt and the high 
Mt slopes) increases with mass as expected from Eq. |2^. 
When the cascade is included, rescattering changes this 
mass dependence. The curvature no longer increases but 
remains approximately constant after the nucleon mass. 



It is useful to compare the flow of the nucleon and the 
fl^ . The nucleon has a smaller mass which, according 
to Eq. decreases the curvature relative to the Q~ . 
However, the nucleon decouples later than the f2^ and 
through hadronic rescattering develops larger transverse 
velocity, which increases the curvature. In the end, the 
fl~ and the nucleon have approximately the same spec- 
tral shape. In contrast, the 4>, which has approximately 
the same mass as the nucleon but which decouples early, 
has very little spectral curvature. To summarize, an in- 
terplay between the differential freezeout dynamics and 
the curved thermal spectra of Eq. results in rich fea- 
tures in the final spectra of different particles. 

Before discussing their impa ct par ameter dependence 
of these rich features (see Sect. IVF ), we study the sen- 
sitivity of the spectra to the EOS. The mass dependence 
of the slopes is a feature of an expanding thermal source 
and differential freezeout. It is not a feature of the under- 
lying EOS. In Fig. ^ the low Mt and high Mt slopes 
are shown for three different EOS. For this discussion, 
the direct comparison of model and data nucleon spectra 
in Fig.p^ may be helpful. For the high Mt slopes for all 

slope 

) the ordering of 



particles, there is a simple ordering, T^^^g > T^^^^l^ > 



T^siopT^ which reflects (through Eq 



the transverse flow, vi^ 



LHS 



> V. 



LHW 



> V. 



LHc 



In the low 



Mt region, the ordering is more complex and reflects 
the space-time structure of the freezeout surface for dif- 
ferent EOS. LHoo evaporates particles shrinking radially 
inward. This causes an enhancement of the particle yield 
at low Mt and gives LHoo a significant slope in the low 
Mt region. Still, with LHS the O— shows much more flow 
at low Mt than it does with LH16 and LHoo indicating 
a large velocity at the end of the mixed phase. 

The curvature in the Mr spectra at small Mt is 
a consequence of the mass dependence of Eq. ^ and 
hadronic rescattering, as discussed above. Now the role 
of hadronic rescattering, or the "pion wind" , is studied 
in detail with Fig. |6[ At low Mt, where the cascade 
is most effective, all the nucleons come from the center 
of the nucleus, as can be seen in Fig.Qb). Accord- 
ingly, the number of collisions is larger and the nucleons 
are accelerated more. At high Mt, the nucleon spec- 



trum (recall Fig. 13) is simply shifted with 2-3 collisions 

by a constant amount, approximately 300 MeV, which 

increases the slope. These collisions happen over a time 

scale of « 10 fm and the collision rate is therefore « — 

5fm 

Collecting these observations, nucleons coming from the 
center of the collision freezeout last, populate the low Mt 
region, and are kicked the most by the pion wind. 



E. Comparison to Central RHIC Data 

Now we compare model predictions to the first RHIC 
spectra. Keep in mind the two m ajor predictions of hy- 



drodynamics introduced in Sect. IV B . First, < Mt > 
should increase significantly. Since LHS was found to give 
the best agreement to SPS flow data, LHS should give 
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FIG. 16: For nucleons in central AuAu collisions, (a) an- 
alyzes the number of collisions experienced by a nucleon as 
function of AIt while (b) analyzes the mean freezeout radius 
of a nucleon as a function of Mt- 



the best agreement at RHIC. Out of all the EOS stud- 
ied, the flow velocity increases the most for LH8. Second, 
the spectra should show the flow profile of Eq. This 
profile is sensitive to the particle mass and flow velocity. 
At RHIC therefore, LH8 predicts a significant change in 
slope from low Mr to high Mt- 

Fig.0 and Fig. 18 show the absolutely normalized 
model spectra for three different EOS compared to data 
for 7T~ ,K~ and p. For each particle type, the transverse 
mass spectrum is strong. Generally, LHoo under-predicts 
the flow profile while LH8 reproduces the spectrum. The 
data indicate a strong macroscopic transverse response, 
as expected of an EOS with speed of sound ~ 1/3. 
Some caveats must be mentioned. It is known that the 
transverse mass spectrum is sensitive to the details of the 
initial profile and longitudinal expansion jl^ . In particu- 
lar, the transverse mass spectrum is modified if the initial 
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FIG. 17: A comparison to PHENIX spectra |72| : (a) com- 
pares 7r~ spectra, (b) compares K~ spectra. Both the model 
and the experimental spectra are absolutely normalized. 



entropy density is distributed according to binary colli- 
sions [Em. However, even when the entropy is distributed 
entirely according to binary collisions, the change in the 
spectrum is small. The strong increase in radial flow from 
the SPS to RHIC is reproduced by the hydrodynamic re- 
sponse of LH8. Thus, prediction (1) is borne out by the 
first spectra at RHIC. Next, look at the shape of the p 
spectrum in Fig^. This flattening at low Mt is charac- 
teristic of a flow profile. We expect a smaller flattening in 
the kaon spectrum. Thus, the rich flow profile of Eq.|25| 
is also borne out in the flrst RHIC data and prediction 
(2) is correct. 

It is worthwhile to plot the tt^ and p spectra on the 
same plot. The spectra almost cross for pT « 2.3 GeV 
It was recently pointed out that the measured 7r~/p ra- 
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FIG. 18: A comparison to PHENIX |[7g| and STAR ^ p 
spectra for (a) low and (b) high Mt, respectively. Both the 
model and the experimental spectra are absolutely normal- 
ized. 



tio is several times above the expected ratio from jet 
fragmentation and from a hydrodynamic calculation that 
does not incorporate chemical freezeout j73|. The ra- 
tio is readily explained in a simple hydro/thermal model 
with additional hadronic scattering. The thermal input 
into RQMD is roughly summarized by Eq. Above 
Mt > 2.0 GeV without rescattering, the slope parame- 
ters of pions and nucleons approach the universal value 

Tsiope ~ 250 MeV. This slope is given by Eq ^ ^ 

the parameters for Tth = 160 MeV and vt = 0.45 c. 



26 with 
Ac- 
counting for hadronic rescattering, the nucleon slope at 
large pt approaches Tsiope ~ 300 MeV and is better de- 
scribed by Tth — 160 MeV and vt — 0.55 c. Hadronic 
rescattering therefore increases the nucleon flow veloc- 
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FIG. 19: A comparison of vr" and p spectra, (a) shows a 
simple thermal model with parameters discussed in the text. 
The spectra are relatively normalized, (b) shows an abso- 
lutely normalized comparison of the complete model spectra 
and PHENIX spectra Q. 



ity slightly, from vt ~ 0.45 c to vt ~ 0.55 c . We then 
adjust hb/T to match the experimental p/p ratio p7| , 
pIp = exp(— 2/XB/r) = 0.65. Then with all the parame- 
ters specified, we draw Eq. |2| for p and tt" in Fig. A 
source expanding with a collective velocity of ~ 0.5 c 
and hadronizing according to a thermal prescription at 
a temperature of T sa 160 MeV, generates the observed 
Tr~ /p ratio once pion nucleon scattering is taken into ac- 
count. 
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The Impact Parameter Dependence of Radial 
Flow 



In peripheral collisions, hydrodynamic features should 
disappear since the mean free path becomes comparable 
to the root mean square radius Rrms- Because ideal hy- 
drodynamics is scale invariant, the hydrodynamic stage 
of the model does not capture finite size effects which 
become increasingly important at larger impact parame- 
ters. However, this does not mean that the radial velocity 
is independent of impact parameter. The hydrodynamic 
lifetime scales approximately as Rrms] the flow velocity 
reflects this lifetime. Finite size effects make the total 
lifetime decrease more quickly than Rrms ■ Finite size ef- 
fects are most important during the freezeout stage which 
is modeled with RQMD and therefore the hydro-|-cascade 
approach can capture some non- trivial features of the im- 
pact parameter dependence. 

In the low Mt region, pion-nucleon scattering in 
RQMD is largely responsible for the w 500 MeV anti- 
proton slope measured by the STAR collaboration [ fzi] . 
In the high Mx region measured by the PHENIX col- 
laboration , the hydrodynamic stage of the model is 
much more significant than hadronic rescattering. Fig. 
shows the low and high Mt slope parameters as a func- 
tion of the number of participants in the collision. For 
each particle, the open symbols show the slope parameter 
without RQMD and the closed symbols show the slope 
parameter with RQMD. For both pious and kaons in the 
low and high Mt regions, only a small impact parameter 
dependence is observed. The pion spectrum is cooled at 
all impact parameters. 

Compare the nucleon and the 4> slope parameters. In 
the low Mt region, shown in Fig.^(a), the nucleon slope 
exhibits a very rapid dependence on the number of partic- 
ipants. Approximately 40% of this slope is due to RQMD 
and the "Hydro Only" curve for the nucleon is relatively 
fiat as a function of impact parameter. The nucleon and 
the (j) have approximately the same mass but the (p is de- 
void of the strong A resonance which drives the flow in 
the nucleon system. Therefore, the "Hydro Only" curve 
for nuclcons is similar to the "Hydro-|-RQMD" curve for 
the (j). 

Next, compare the nucleon and the (j) slopes in the high 
Mt region, shown in Fig.|2^(b). Here ttN — + A is less 
significant and RQMD is responsible for only « 15% of 
the nucleon slope. Consequently, RQMD makes up the 
small dilference between the nucleon and the (f> "Hydro 
Only" curves, and the final slope parameters of the two 
particles are similar in shape and magnitude. Experi- 
mentally, the difference in slope parameters between the 
nucleons and the ((> can be used to assess the contribution 
of the hadronic phase in the low Mt region. 

Turning to the experimental data, we first examine the 
low Mt region as measured by the STAR collaboration. 
We compare the model b-dependence of the anti-proton 
and K~ slope parameters to the experimental data in 
Fig. E| and E^. Comparing the anti-proton slopes for 
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FIG. 20: The (a) low Mr and (b) high Mr slope parameters 
(see Fig. ^ as a function of the number of participants in 
the collision. 



different EOS, we see that the rapid b-dependence is not 
a consequence of a change in EOS. Note that LHoo has 
a larger slope at small Mt than LH8. This is an artifact 
of the exponential fit. The spectrum for LHcx) in central 
collisions is shown in Fig|l^ and is not exponential. The 
measured spectra are well described by an exponential 
in this region |Q. Coarsely, the model reproduces the 
strong b-dependence of the anti-protons and the weaker 
b-dependence of the kaons. However, for more periph- 
eral collisions, the measured slope parameters fall some- 
what faster than the model predicts. Naively, this indi- 
cates that in the most peripheral collisions the hydrody- 
namic description is only beginning to work and finite 
size effects (e.g viscosity) should be taken into account 
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FIG. 21: Model and preliminary STAR K' jTll slope pa- 
rameters in the fit range < Mt — M < 0.5 GeV as a function 
of the number of participants for different EOS. (a) compares 
model predictions for each EOS with and without the RQMD 
after-burner; (b) compares the model to data with the RQMD 
after-burner. 



FIG. 22: Model and prehminary STAR p Q slope param- 
eters in the fit range < Mt — M < 0.5 GeV as a function 
of the number of participants for different EOS. (a) compares 
model predictions for each EOS with and without the RQMD 
after-burner; (b) compares the model to data with the RQMD 
after-burner. 



in the QGP and mixed phases which are modeled with 
the hydrodynamics. The "Hydro Only" curves presented 
in Fig. ^ and ^ definitely do not reproduce the strong 
b-dependence of the slopes. With RQMD on the other 
hand, finite size effects during the late hadronic stages 
are modeled and most of the rapid b-dependence seen in 
the data is reproduced. 

Turning to the high Mt region, the absolutely nor- 
malized model spectra are compared directly to the ab- 
solutely normalized PHENIX spectra for different cen- 
tralities and particles in Fig.|2^. For and p the spec- 



trum is reproduced in the most central bin (0-5% central 
with < Np >— 350) and in the semi-peripheral bin (30- 
60% central with < Np >~ 76). For pions, the spectral 
shape in the most central bin is reproduced. However, in 
the semi-peripheral bin, the pion spectrum resembles a 
power-law rather than a thermal spectrum. This change 
in shape from peripheral to central has been attributed 
to jet-quenching [^, [75| . 
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FIG. 23: The absolutely normalized model spectra for LH8 
compared to PHENIX spectra |72| at different centralities. 
(a) is for vr^, (b) is for and, (c) is for p. 



V. ELLIPTIC FLOW FROM THE SPS TO RHIC 

In non-central collisions the particles emerge with an 
elliptic flow. The spectator matter flies down the beam 
pipe and the excited nuclear matter is formed in the 
transverse plane with an almond shaped distribution. 
Subsequently, if pressure develops in the system, the pres- 
sure gradients are larger in the impact parameter direc- 
tion (the x-direction) than in the y-direction. Then, the 
excited matter expands preferentially in the x-direction. 
The magnitude of this elliptic response is quantified ex- 
perimentally by expanding the distributions in a Fourier 
series 



dN 



dN 



Pt dpT dy I 



2tt pTdprdy 



+...) 



{l + 2v2{pT,y) cos{2(j)) 
(27) 



where is measured around the z-axis relative to the 
impact parameter, which points in the x direction. The 
elliptic flow, V2{pT,y) = (cos(2(/)))py ^y, gives a measure 
of the dynamic response of the excited nuclear matter to 
the initial anisotropy. 

The initial spatial anisotropy is quantified using the 
parameter e (see Eq. P). The hydrodynamic response 
is linear in e [ p9[ and therefore V2 is sometimes divided 
by e to compare different impact parameters and nuclei 
p2| , [76[ . As the system expands, the eccentricity e de- 
creases. Since e is the driving force behind the elliptic 
flow, the elliptic development finishes before the radial 
development. Therefore, elliptic flow is generated by the 
early pressure, although this statement must be qualified 
(see below). The spatial anisotropy that remains after 
the colhsion is quantified by S2 (see Eq. |o|) . S2 measures 
how much of the initial spatial anisotropy e was not used 
during the collision for the production of elliptic flow. 



A. Qualitative Changes from the SPS to RHIC 



In Sect. IV A radial flow was used to constrain the 



EOS. The best (though certainly not unique) description 
of the data was given by LH8. For consistency, we require 
the same EOS to describe the elliptic flow data. 

Fig.^(a) shows integrated elliptic flow of pions as 
a function of the total multiplicity for different EOS. 
Fig. |2^(b) shows the relative contribution of RQMD to 
the integrated elliptic flow. Note that elliptic flow in- 
creases for all EOS and dramatically so for LH8. As- 
sume momentarily that elliptic flow for Hydro-|-RQMD 
stops developing at a temperature of T « Tc « 165 MeV. 
(Note however that the radial flow develops well be- 
low this temperature). The dramatic increase of ellip- 
tic flow in Fig.^(a) can be understood as the dynamic 
response of the QGP pressure. Recall Fig. ||(b), which 
plots anisotropy of the hydrodynamic stress tensor versus 
time and pay particular attention to the T « Tc points 
(the solid symbols) on the LH8 and RG curves (ignore 



22 




0.02 



0.01 



^^LH8 EoS 
^ LHI6E0S 
- LHco EoS 
* Hydro T,=1 20 MeV 
o RG EoS 
' Na49+STAR Data 



0.05 




Hydro+RQMD 
Hydro Only 

=120 MeV 
Hydro+RQMD 



0.01 



0.05 



LH8 Hydro+RQMD 
LH8 Hydro Only 
Hydro T,=1 20 MeV 
-LHoo Hydro+RQMD 
-LH»= Hydro Only 



CO 
c 
o 

Q. 



-0.05 




200 \ 
SPS 



406\ 
RHIC 



600 800 
dN,,/dy 



FIG. 24: Panels (a)-(c) show three related quantities as 
a function of the total multiplicity in a PbPb collision at 
b = 6fm. (a) shows the integrated elliptic flow V2 of pi- 
ous for different EOS and freezeout conditions; (b) shows the 
integrated elliptic flow V2 with and without hadronic rescat- 
tering; (c) shows the spatial anisotropy S2 with and without 
hadronic rescattering. At the SPS, the NA49 112 data point 
is extrapolated to b=6fm using Fig. 3 in |^^. At RHIC, the 
STAR V2 data point is extrapolated to Nch/NJ^''^ = 0.545 
(b=6fm in AuAu) using Fig. 3 in jTSl]. 



LHoo for now). At the SPS, the anisotropy of the stress 
tensor increases rapidly at first and then stalls. The final 
stress tensor anisotropy is small at the end of the mixed 
phase. At RHIC, elliptic flow develops more rapidly and 
stalls only when the anisotropy is large. Thus, provided 
the elliptic flow stops developing at Tc, the elliptic flow 
increases dramatically as the QGP pressure appears. For 
a KG EOS at the SPS, there is no mixed phase and no 
stall and consequently the RG elliptic flow is signiflcantly 
larger than LH8 and the data. However with RHIC col- 
lision energies, LH8 begins to behave as an ideal QGP. 
Consequently, at RHIC the RG elliptic flow is only 20% 
larger than that of LH8 and of the data. 

To understand when elliptic flow stops developing, it 
is important to track the spatial geometry of the under- 
lying source. When the spatial anisotropy, S2, is neg- 
ative, the pressure drives elliptic flow. However, as S2 
approaches zero, the pressure gradients drive radial mo- 
tion rather than elliptic motion. The elliptic development 
then stops. Fig. p^(c) shows the spatial anisotropy, S2, for 
pions as a function of multiplicity from the SPS to RHIC. 
Compare the LH8 curves (the stars, the solid squares, 
and the open squares) seen in Fig.|2j(c). The open 
squares (LH8 Hydro Only) show the spatial anisotropy, 
S2, at the end of the mixed phase, the closed squares 
(LH8 Hydro-|-RQMD) show S2 after the cascade, and the 
stars show S2 when the hydrodynamic evolution is contin- 
ued to Tf 120 MeV. After the mixed phase (LH8 Hy- 
dro Only), the matter retains some of its initial almond 
shape. Continuing the hydrodynamics destroys the ini- 
tial almond shape completely and increases V2 by a factor 
of « 2 — 3 (see Fig. |^(c)). Cascading also changes the 
almond shape but increases V2 by only a factor of « 1.5. 
In cither case, S2 crosses zero between the SPS and RHIC 
and therefore elliptic development in the hadronic stage 
ceases to be significant between the SPS and RHIC. 

This fact is illustrated with Fig. ^ which contrasts the 
RQMD contribution to V2{pt) with the contribution of 
the hadronic phase in a pure hydrodynamic calculation 
at the SPS and RHIC. At the SPS, V2{pt) increases by 
a factor of two when the the hydrodynamics is continued 
to Tf — 120 MeV. When the hydrodynamics is replaced 
with RQMD, V2{pt) also develops but only by approx- 
imately 20%. At RHIC, the spatial asymmetry is com- 
pletely destroyed by the end of the mixed phase and the 
elliptic development is frozen for all pt ■ Continuing with 
the cascade or the hydrodynamics increases the elliptic 
flow marginally. 



B. The Impact Parameter Dependence of Elliptic 
Flow 

Fig.^ shows the b or Np dependence of integrated 
pion elliptic flow at the SPS and RHIC. 

The data restrict the underlying EOS. At the SPS, the 
data favor a soft EOS - LH > 0.8GeV/fm^ LH4 and 
RG EOS generate too much elliptic flow. For LH8-LH16, 
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FIG. 25: The dependence of V2{pt) on the switching tem- 
perature Tswitch at (a) the SPS and (b) RHIC for an impact 
parameter of b = 6 fm. 

the model is about 20% above the data. However, the 
model to data comparison is not completely fair - the 
data points are integrated over rapidity, while the model 
points are only strictly valid for mid-rapidity. This prob- 
ably accounts for the residual model/data discrepancy. 
As the latent heat is increased beyond LH32 to LHoo , 
the elliptic flow begins to rise. The origin of this ellip- 
tic flow was described in and results from the slow 
evaporation of particles in an asymmetrical fashion over 
a long time. This elliptic flow is generated without ra- 
dial flow and the px dependence of V2 for nucleons 
(see below)is modified accordingly At RHIC (Fig. 
p6|(b)), the comparison is fair and the data again favor 
a relatively soft EOS, LH8-LH16. Thus the elliptic flow 



FIG. 26: (a) V2 for pions at the SPS as a function of par- 
ticipants (relative to the maximum) compared to NA49 data 
[{ttI . (b) V2 for charged particles at RHIC as a function of 
participants (relative to the maximum) compared to STAR 
data [ [r8| . 

data at the SPS and RHIC are consistent with a single 
underlying EOS. 

Note that the ordering of the EOS in Fig^ differs at 
the SPS and RHIC. At the SPS, LH8 and LH16 gen- 
erate approximately the same elliptic flow. At RHIC, 
the hard QGP phase lives for substantially longer with 
LH8 than with LH16 and therefore generates more ellip- 
tic flow. Additionally at RHIC, LH4 generates more ellip- 
tic flow than a RG EOS. Thus, the elliptic flow indicates 
that at high energy densities LH4 (with w 1/3) has a 
larger speed of sound than a RG EOS (with ~ 1/5). 
At asymptotically, high energy densities all EOS in the 
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LH(x) family approach the massless ideal gas limit. 



C. The pt Dependence of Elliptic Flow 

Having discussed qualitative changes from the SPS to 
RHIC, we explore the pt dependence of elliptic flow. Ex- 
perimental measurements are performed over a range of 
impact parameters. To find V2{pt) in a specific impact 
parameter range, 6mm < b < bmax, the following inte- 
grals must be performed, 



K:::MPT,y;b)^{b)2nbdb 



dN 
dy dpq 



■(&) 2TTbdb 



(28) 



Again, we drop the y, bmin and bmax labels below when 
it is not confusing. (pt)™*""'""'* denotes the elliptic 
flow integrated over all events, or V2{pt)o^ ■ 

Fig. |^(a), (b) and (c) show ^;2(pt)"""''"'" for nega- 
tive hadrons, kaons, and nucleons at RHIC. Look first at 
the negative hadrons (a): Although LH8 and LH16 both 
show a strong linear rise, the slope is smaller for LH16. 
For LHoo, V2{pt) is curved and bends over. For small p^, 
LHcxD is above LH8, but by px ~ 2.0 GeV, LHcx) is sub- 
stantially below LH8. The data show a strong linear rise 
and agree remarkably well with the slope of LH8. V2{pt) 
slightly favors LH8 over LH16. The kaon V2{pt) curve 
has the same shape and magnitude as the h~ spectrum. 
At the SPS the kaon elliptic flow is slightly negative 
This anti-elliptic flow is most likely a remnant of the re- 
pulsive mean field observed at higher baryon densities. 
At RHIC, the baryon density is lower than at the SPS 
and kaons should flow along with the pions if the space 
time picture of the model is correct. 

For nucleons, the V2 spectral shape is different and is 
initially curved upwards. A useful thermal model has 
been given to explain the shape of V2{pt) 0. For 
nucleons, LH8 and LH16 are concave up, indicating a 
strong radial expansion. By contrast, LHoo shows a lin- 
ear rise in V2{pt)i indicating a weak transverse expan- 
sion. As discussed in Sect. |l^ LHoo slowly evaporates 
particles into RQMD and generates elliptic flow only at 
small Pt- The curvature of V2 (pt) for LHoo resembles the 
Pt dependence expected if only surface evaporation were 
present However, LHoo does develop a substantial 
radial flow over its long lifetime which gives the LHoo 
V2{pt) some shape. The data favor the strong transverse 
expansion of LH8 over the weak expansion of LHoo . 

We now demonstrate that pion nucleon scattering on 
top of a baseline elliptic flow is responsible for the cur- 
vature of V2{pt) seen in data. Fig. p8| shows V2{pt) with 
and without hadronic rescattering. Here the discussion 
parallels the discussion of the radial flow. Pion nucleon 
scattering increases the radial flow of the nucleon spec- 
trum and cools the pion spectrum. Consequently the 
pion V2 spectrum with the RQMD after-burner is slightly 
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FIG. 27: i;2(pt)'"'""""" (see Eq. g) for three different EOS 
compared to data. Panels (a), (b) and (c) are for negative 
hadrons, all kaons, and p-l-p, respectively. The data are from 
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FIG. 28: 112 as a function of pr for (a) all pions and (b) 
nucleons+hyperons with and without the RQMD after-burner 
for AuAu collisions at b=6fm, for LH8 EOS. 



FIG. 29: V2{pt)/^ for semi-central (b=4fm) and semi- 
peripheral (b=10fm) collisions with LH8 EOS. Panel (a) is 
for pions and (b) is for nucleons-|-hyperons. 



above the "Hydro Only" spectrum. Similarly, the nu- 
cleon V2 spectrum with the after-burner is curved upward 
by TT — iV scattering within RQMD. Similar features were 
found in all the EOS studied above. 

Now to illustrate the impact parameter dependence, 
Fig. ^ shows the b-dependence of V2{pt)- To compare 
different impact parameters, V2{pt) for pions and nucle- 
ons is divided by the initial space anisotropy e, for periph- 
eral (b=10fm) and semi-central (b=4fm) collisions. The 
model response basically follows naive geometric consid- 
erations. However, closer inspection reveals that the 
model captures some finite size effects during the late 
hadronic stages. As the impact parameter is scanned, 
the total number of pions decreases roughly as the num- 



ber of participants, and pion-nucleon scattering decreases 
similarly. Consequently, for central collisions pions show 
slightly larger elliptic flow at small pT while nucleons 
show a smaller (more curved) elliptic flow at small pt- 
Thus, together these curves indicate a slightly stronger 
hadronic expansion in central collisions. 

We now return to the SPS and compare the model 
to NA49 elliptic flow data. Fig. ^ compares the model 
and data V2{pt) for pions and nucleons. These data are 
generally less well produced than at RHIC. Some caveats 
should be mentioned. The data are forward in rapidity, 
4 < ?/ < 5, while strictly speaking, the model is only valid 
for mid-rapidity data (y=2.92). For pions, V2 has a fairly 
flat rapidity dependence but the total abundance changes 
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FIG. 30; 112 as a function of pr for pions and protons in 
PbPb collisions at the SPS. The data points are for events 
with b ~ 6.5 — 8fm or more precisely for events with 45- 
60% of the beam energy in the NA49 zero-degree calorimeter 
[^ . The model points are the results of integrals between 
b = 6.5 — 8 fm or V2{pT)t'.5 as described by Eq. ESl The model 
data comparison is not completely fair - the model points are 
for mid-rapidity while the data points have been integrated 
over rapidity (see text). 



significantly from y=2.92 to y=5. For nucleons, V2 has 
strong rapidity dependence and is almost a factor of two 
larger at mid-rapidity. At mid-rapidity, the elliptic flow is 
certainly stronger, which should improve the agreement 
with the model. Furthermore, in the forward rapidity 
region, the dynamics are complex; pions and nucleons 
have a significant directed flow indicating that details 
of stopping may play a significant role. In conclusion. 



much more data are needed to establish the applicability 
of hydrodynamics at the SPS. 



D. Summary and Comparison 

The principle motivation of this work was to demon- 
strate that the body of heavy ion data from the SPS to 
RHIC can be described with thermodynamics and hy- 
drodynamics. To this end, we have compared of the hy- 
dro-|-cascade model of |^ to the radial and elliptic flow 
data from the SPS and RHIC. A simultaneous analysis 
of available flow constrains the EOS. Only an EOS ex- 
hibiting the hard and soft features of the QCD phase 
transition systematically reproduces the observed radial 
and elliptic flow. 

The model incorporates strong radial and elliptic flow, 
chemical freezeout at the phase boundary, subsequent 
hadronic rescattering and differential freezeout. With 
these ingredients, the model explains a number of fea- 
tures of the new RHIC data. 

• The "anomalous" p/tt~ ratio (which exceeds one 
for pt ^ 2 GeV) is a simple consequence of the in- 
crease in the radial flow and chemical freezeout. In 
a thermal picture, anti-protons are enhanced rela- 
tive to proton-proton collisions. Then, the strong 
radial flow drives these anti-protons to large px- 
Subsequent hadronic rescattering makes the spec- 
tra cross. 

• The Mt spectra (which are "curved" as opposed 
to simply exponential) are readily explained in a 
hydro-|-cascade model. The curvature is due to 
a combination of the flow profile expected from 
hydrodynamics |5^, ^ and hadronic rescattering. 
With these ingredients, the mass dependence of the 
spectra measured by the STAR and PHENIX col- 
laborations are naturally explained. The strong b- 
dependence of the STAR slope parameters for anti- 
protons is a consequence pion-nucleon scatter- 
ing. In contrast, the slope parameter for the 0, 
which suffers few hadronic collisions in the model, 
shows little b-dependence. 

• The observed elliptic flow rises rapidly as a function 
of pt and favors a strong transverse expansion. Un- 
like the radial spectra, the elliptic spectra are less 
sensitive to hadronic rescattering and differential 
freezeout. Therefore, our results on the px spec- 
trum of V2 are largely similar to the hydrodynamic 
analysis in pO| . 

These features are generic to a radially and elliptically 
expanding thermal source and do not immediately impli- 
cate hydrodynamics as the dynamic origin of the radial 
and elliptic flow. However, running hydrodynamic up to 
the phase boundary (with the same EOS) quantitatively 
reproduces the necessary radial and elliptic flow velocities 
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both at the SPS and RHIC. In particular, hydrodynamics 
reproduces the observed changes from the SPS to RHIC: 

• In a hydrodynamic framework, the radial flow ve- 
locity increases at high energy densities for an EOS 
with a phase transition to the QGP |2^. At 
the SPS, LH8 gave the best agreement with spec- 
tra and predicted a « 20% increase in the radial 
flow velocity from the SPS to RHIC |l|. The 
flrst RHIC spectra conflrm the predictions of LH8 
and the predictions of other hydrodynamic works 
[ p^ , [20[ |. Generally, LH8 has a smaller latent heat 
than that used in other works and therefore LH8 
predicts a larger increase in the radial flow. In par- 
ticular, already at RHIC, the rriT spectrum of the 
r2~ is significantly curved by the radial flow |32 . 



• At RHIC and to a lesser extent at the SPS, the mag- 
nitude of the integrated elliptic flow is reproduced. 
During the early stage of this work, a ~ 40% change 
in elliptic flow the SPS to RHIC was predicted and 
subsequently observed |7|]. This increase is a di- 
rect consequence of the QGP pressure [|l^, ^ and 
the early freezeout of elliptic flow at the SPS pl|. 



E. Fixing the EOS 



F. Outlook 

Finally, we turn to open problems. Hanbury-Brown 
Twiss (HBT) correlations provide information about the 
spatial and temporal extent of the freezeout region. Such 
measurements have been performed at the AGS, the SPS 
lH and recently at RHIC Q. Although the HBT radii 
fall with the pion pair momentum Kt providing addi- 
tional evidence for transverse flow |Q, the measured 
radii at RHIC are approximately only 50% percent of 
our preliminary radii |85|] . 

The dynamic origin of these small HBT radii is not un- 
derstood and much more work is needed |]86| . The small 
radii indicate that although the final velocities are cor- 
rectly reproduced within the model, the model expansion 
time is too long. Future experiments will measure HBT 
radii and deuteron production as a function of impact pa- 
rameter and azimuthal angle. Such detailed experimental 
information is essential to a complete understanding of 
the excited matter produced in ultra-relativistic heavy 
ion reactions. 
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Taking the radial and elliptic flow together, we argue 
that the momentum correlations in the data reflect the 
hydrodynamic response of excited matter exhibiting the 
soft and hard features of the QCD phase transition. For 
an EOS without softness, e.g. a resonance gas EOS, the 
flow of multi-strange baryons is dramatically missed (see 
Fig. [T^ ). In addition, the elliptic flow is significantly too 
large both at the SPS and at RHIC (see Fig. ^(a) and 
(b)). These observations indicate that without softness 
the initial hydrodynamic response of the fireball is too 
strong. 

For an EOS without a hard component, e.g. LHcxd, 
the spectra are significantly too soft both at the SPS 
and RHIC (see Fig. |l] and Fig. |l|). The flow of the 
multi-strange baryons is even too small (see Fig. |l2|) . Al- 
though LHoo generates a large V2 by evaporating parti- 
cles anisotropically through the freezeout surface, the pt 
dependence of this elliptic flow is qualitatively wrong (see 
Fig. ^7|). For LHoo, hadronic rescattering does generate 
some transverse flow, but this transverse flow is insuffi- 
cient to explain the strong px dependence of the elliptic 
fiow. The strong curvature for V2{pt) seen in the nucleon 
data implicates a violent explosion and this violent ex- 
plosion is not generated by LHoo. Out of all the EOS 
considered, the best agreement is found with LH8. The 
same EOS was deduced prior to the analysis of the first 
RHIC data (|l). LH8 has a latent heat of 800 MeV/fm^ 



and represents a balance between soft and hard. 



[4; 

[5 
[6 

[r 

[8 
[9 

[lo; 
[11 

[12 

[is: 

[14 
[15 

[16" 



REFERENCES 

Quark Matter 1997, Nucl. Phys. A638, (1998). 
Quark Matter 1999, Nucl. Phys. A661(1999). 
Quark Matter 2001, to be published in Nucl. Phys. 
A, Stony Brook, NY, January 15-20 (2001) ; 
http://www.rhic.bnl.gov/qm2001/program.html 

E. V. Shuryak, Phys. Reports C61, 71 (1980). 

e.g., M. Oevers, F. Karsch, E. Laermann, and P. 
Schmidt, Nucl. Phys. Proc. Suppl. 73, 465 (1999). 
P. Braun-Munziger, J. Stachel, J. P. Wessels, N. Xu, 
Phys. Lett B344, 43 (1995). 

F. Becattini, M. Gadzdzicki and J. SoUfrank, Eur. 
Phys. J. C5, 143 (1998) 

J. SoUfrank, J. Phys. G23, 1903 (1997). 
J. Cleymans and K. Redlich, Phys. Rev. C59, 354 
(1999) ; J. Cleymans, H. Oeschler and K. Redlich, 
Phys. Rev.C59, 1663 (1999) ; 

E.V. Shuryak and O.V. Zhirov, Phys. Lett. B89, 253 
(1979). 

e.g., R. Stock, in Quark Matter 1999, Nucl. Phys. 
A661, 419c (1999). 

e.g., U. Heinz, Nucl. Phys. A685, 414 (2001) and 
references therein. 

M. Gyulassy and X.N. Wang, Comp. Phys. Comm. 
83, 307 (1994); Phys. Rev. D44, 3501 (1991). 
D. Molnar and M. Gyulassy, nucl-th/0104073. 
M. Bleicher and H. Stocker, hep-ph/0006147. 
Josef SoUfrank et al, Phys. Rev C55, 392 (1997). 



28 



U. Ornik, M. Pluemer, B.R. Schlei, D. Strottman, [50 
R. M. Weiner Phys. Rev. C54 1381 (1996). 
P.P. Kolb, J. SoUfrank, U. Heinz, Phys. Rev. C62, [51 
054909 (2000). 

P.P. Kolb, P.Huovinen, U. Heinz, H. Heiselberg, [52 
Phys. Lett B500, 232 (2001). 

P. Huovinen, P.P. Kolb, U. Heinz, H. Heiselberg, [53 
Phys. Lett B503, 58 (2001). 

D. Teaney, J. Lauret, and E.V. Shuryak, Phys. Rev. [54 
Lett. 86, 4783 (2001). 

CM. Hung, E.V. Shuryak, Phys. Rev. Lett. 75, 4003 [55 

(1995) . 

C. Greiner, D. H. Rischkc, Phys. Rev. C54 1360 [56 

(1996) . 

NA49 Collaboration, T. Alber et ai, Phys. Rev. [57 
Lett. 75, 3814 (1995). [58 
PHOBOS Collaboration, B.B. Back et al, Phys. Rev. 
Lett. 85, 3100 (2000) ; [59 
L. Van Hove, Z. Phys. C21, 93 (1983). 
J. Blaizot and J.-Y. OUitrault, Phys. Rev. D36, 916 [60 
(1987); [61 
M. Kataja et ai, Phys. Rev. D34, 794 and 2755 [62 
(1986). [63 
J.-Y. OUitrault, Phys. Rev. D46, 229 (1992) [64 

D. H. Rischke et ai. Heavy Ion Physics 1, 309 (1995). 
nucl-th/9505014. [65 
U. Heinz, Nucl. Phys. A661, 141c (1999). 

S. Bass and A. Dumitru, Phys. Rev. C61, 064909 [66 
(2000). 

M. Belacem et ai, Phys. Rev. C58, 1727 (1998). [67 
H. van Hecke, H. Sorge, N. Xu, Phys. Rev. Lett. 81, 
5764 (1998). [68 
P.P. Kolb, J. SoUfrank, U. Heinz, Phys. Lett. B459, 
667 (1999). 

H. Sorge, Phys. Rev. C 52, 3291 (1995). 

J.D. Bjorken, Phys. Rev. D27, 140 (1983). [69 

D. Teaney, in preparation. 

F. Cooper and G. Frye, Phys. Rev. DIO, 186 (1974). [70; 
L.D. Landau and E.M. Lifschitz, Fluid Mechanics, 
2nd Edition (Butterworth and Heinemann,1987) [71 
M. Gyulassy and T. Matsui, Phys. Rev. D29, 419 
(1984). 

R.J. LeVeque, Numerical Methods for Conservation 
Laws, (Birkhauser-Verlag,1990). [72 
V. Schneider et al, J. Comput. Phys. 105, 92 (1993). 
Dirk H. Rischke, Yaris Pusun, Joachim A. Maruhn, [73 
Nucl. Phys. A595, 346 (1995) ; Nucl. Phys. A595, [74^ 
383 (1995). [75 
P.P. Kolb, U. Heinz, P. Huovinen, K.J. Eskola, K. 
Tuominen, hep-ph/0103234. [76 
P. Jacobs and G.Cooper, STAR SN402(1999). 
STAR CoUaboration, C. Adler et ai, Phys. Rev. Lett [77 
86, 4778 (2001). 

PHOBOS Collaboration, B.B. Back et al, submitted [78 
to Phys. Rev. Lett., 

http://www.phobos.bnl.gov/Publications/Ph_pub.htm [79 

PHENIX CoUaboration, K. Adcox et al, Phys. Rev. 

Lett. 86, 3500 (2001). [80 



PHOBOS Collaboration, B.B. Back et al, submitted 
to Phys. Rev. Lett, nucl-ex/0105011. 
S. Soff, S.A. Bass, and A. Dumitru, Phys. Rev. Lett 
86, 3981 (2001). 

R. Venugopalan and M. Prakash, Nucl. Phys. A546, 
718 (1992). 

J.P. Blaizot and J.-Y. OUitrauU, Nucl. Phys. A458, 
745 (1986). 

S.M. Johns, P.J. Ellis, J.M. Lattimer, Astrophys. J. 
473, 1020 (1996). 

F. Cooper, G. Frye, E. Schonberg, Phys. Rev. Dll, 
192 (1975). 

Cs. Anderhk et ai, Phys. Rev. C59, 388 (1999) and 
references therein. 

K.A. Bugaev, Nucl. Phys. A606, 559 (1996). 

D. Teaney, Ph.D. Thesis, State University of New 
York at Stony Brook (2001). 

E. Schnederman, J. SoUfrank, U. Heinz, Phys. Rev. 
C48, 2462 (1993). 

H. Sorge, Phys. Rev. Lett. 78, 2309 (1997). 
H. Sorge, Phys. Lett. B402, 251 (1997). 
H. Sorge, Phys. Rev. Lett 82, 2048 (1999). 

H. Sorge, Phys. Rev. Lett. 82, 2048 (1999). 
STAR Collaboration, C. Adler et al. , nucl- 
ex/0107003. Submitted to Phys. Rev. Lett. 

CM. Hung, and E. Shuryak, Phys. Rev. C57, 1891 
(1998). 

NA49 Collaboration, H. Appelshauser et ai, Phys. 
Rev. Lett. 82, 2471 (1999). 

D. H. Rischke and M. Gyulassy, Nucl. Phys. A608, 
479 (1996). 

E. Andersen et al., WA97 collaboration, Phys. Lett. 
B433, 209 (1998); R. Lietava for the WA97 collab- 
oration in Proc. Strangeness 98, J. Phys. G: Nucl. 
Part. Phys. 25, (1999). 

G. Roland for the NA49 CoUaboration, Nucl. Phys. 
A638, 91c (1999). 

P.J. Siemens and J.O. Rasmussen, Phys. Rev. Lett. 
42, 880 (1979). 

M. Calderon, for the STAR Collaboration, at Quark 
Matter 2001 [§; STAR CoUaboration, C. Adler et 
al, to be submitted to Phys. Rev. Lett. 
http://www.star.bnl.gov/ . 

J. Velkovska, for the PHENIX Collaboration, at 
Quark Matter 2001 [§; nucl-ex/0105012. 

I. Vitev and M. Gyulassy, nucl-th/0 104066. 
X. Wang, at Quark Matter 2001 |3|. 

G. David, for the PHENIX Collaboration, at Quark 
Matter 2001 @ 

S.A. Voloshin and A.M. Poskanzer, Phys. Lett. 
B474, 27 (2000). 

A.M. Poskanzer and S.A. Voloshin for the NA49 Col- 
laboration, Nucl. Phys. A661, 341c (1999). 
STAR Collaboration, K.H. Ackermann et ai, Phys. 
Rev. Lett 86, 402 (2001). 

P. Huovinen, P.P. Kolb, U. Heinz, in Quark Matter 
2001; nucl-th/0104020 

Shunji Nishimura for the WA98 Collaboration, Nucl. 



29 



Phys. A661, 464c (1999). 
[81] NA49 Collaboration, H. Appclshauser et ai, Phys. 

Rev. Lett. 80, 4136 (1998). 
[82] E.V. Shuryak in QM '99 , Nucl. Phys. A661, 119c 

(1999); 

[83] e.g., U. Heinz and B.V. Jacak, Ann. Rev. Nucl. Part. 
Sci, 49, (1999). 



[84] STAR Collaboration, C. Adler et ai, Phys. Rev. 

Lett. 87, 082301 (2001). 
[85] D. Teaney, RHIC-INT Workshop 2001. Lawrence 

Berkley National Laboratory, May 2001. 

http://www-rnc.lbl.gov/ nxu/workshop/rwwOl.html 
[86] e.g., Tetsufumi Hirano et al, nucl-th/01100G9. 



